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Abstract 

The ratio of J/ip over Drell-Yan cross-sections as a function of the energy of 
the zero degree calorimeter in Pb Pb collisions has been computed in a comovers 
interaction model and compared with the results of a deconfining model. The predic- 
tions of the two models are different both for peripheral events and for very central 
ones. These differences are analyzed and the results of the models confronted with 
available data - not only for Pb Pb but also for pA and SU collisions. 



1 Introduction 

The NA50 interpretation of the data on J / ip suppression in pA, SU and Pb Pb collisions 
is as follows I]. The pA, SU and peripheral Pb Pb data (up to Et ~ 35 -j- 40 GeV), 
are well described with nuclear absorption alone with an absorptive cross-section a a b s = 
6.4 ± 0.8 mb. At Et ~ 40 GeV there is a sudden onset of anomalous suppression. A 
second accident occurs at Et ~ 90 -j- 100 GeV, close to the knee of the Et distribution, 
where a change of curvature in the shape of the suppression is observed - followed by a 
steep fall-off. However, at variance with this view, the most peripheral point in Pb Pb has 
consistently lied above the NA50 nuclear absorption curve - which extrapolates pA and 
SU data. This tendency is now confirmed by the preliminary data on the J/ip suppression 
versus the energy, Ezdc, of the zero degree calorimeter. Here several peripheral points 
lie above the NA50 nuclear absorption curve and, more important, they exhibit a steeper 
suppression pattern. It is, therefore, important to study these data and to compare them 
with the results of deconfining and comover interactions models[|. This comparison is also 
of interest in the large Et region. Indeed, the NA50 data for Pb Pb collisions at large Et 
can be described either in a deconfining model or as a result of the interaction with 
comovers 0. 

Our purpose here is to confront with each other the results of two specific models, 
namely a deconfining and a comovers interaction || model - and to compare them 
with experiment. In the first model |4j], one requires either two sharp thresholds, as in 
|J, or a single threshold smoothed with an arbitrary function. Both in this model and 
in the comovers one 0, it is necessary to introduce the fluctuations in Et 01 m order to 
reproduce the fall-off of the data at large E T (beyond the "knee" of the E T distribution). 
It has been shown in 0] that the effect of the Et fluctuations is significantly larger in 
the deconfining scenario than in the comovers approach. This leads to agreement with 
the large Et data in the first case - while the fall-off obtained in the comovers approach 
is too weak || [7|. However, all data beyond the knee have been obtained in the so- 
called Minimum Bias {MB) analysis, where only the ratio of J/ip to MB cross-sections 
is measured and it is divided by a theoretical ratio of DY to MB cross-sections, i.e. 



In the theoretical model used by the NA50 collaboration ||, the ratio (DY/MB) t h satu- 
*For reviews of deconfining and comover interaction models, see 0. Alternative models have also 
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been proposed ||. 
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rates at large Et- In contrast, it has been argued in || that this ratio also falls at large 
Et- This is due to the shift in Et between the DY and MB event samples, resulting 
from the Et taken by the dimuon. This shift is very small (of the order of a few per 
mil) and, therefore, has no visible effect up to the knee of the Et distribution - where 
this distribution is rather flat. However, it does have a sizeable effect in the tail, due 
to the very steep fall-off of the distribution. Actually, experimental data on the ratio 
DY/MB are available ||. They do show a fall-off at large Et-, but the statistical errors 
are large. Obviously, the same effect should be present in the ratio J/ip over MB and 
would cancel out in the true ratio of J/ip to DY cross-sections - i.e. the one obtained in 
the so-called standard analysis. However, it does not cancel in Rmb because this effect is 
not included in the theoretical model for the last factor in the r.h.s. of Eq. ([!]), used by 
the NA50 collaboration |§. When this effect, as estimated in ||, is taken into account, 
the comovers approach does describe the data for Rmb-, whereas, in the deconfining model 
[|]], the fall-off beyond the knee is too strong. This is shown in Fig. 1 where the dashed 
lines are the results in the comovers model || with and without Et loss. The upper solid 
line is the result obtained in the deconfining model [|J. In this model one calculates the 
true ratio of J/ip to DY differential cross-sections - not Rmb- Therefore, this line should 
be compared with the results in the standard analysis. To compare with the data in Fig. 
1, which where obtained with the MB analysis, one should take into account the Et loss. 
Incorporating this effect, following the prescription ||, one obtains the lower full curve in 
Fig. 1. 

Clearly, a good way to avoid this problem is to measure Rmb as a function of E Z dc- 
Indeed, the main contribution to Ezdc is due to the energy of the spectator nucleons - 
which is not affected by the presence or absence of the dimuon trigger. 

The plan of this paper is as follows. In Section 2 we derive the Et - Ezdc correlation. 
It turns out that this correlation allows to relate to each other the MB experimental 
distributions in the two variables - Et and Ezdc- In Section 3 we compute the J/ip 
suppression versus Ezdc and compare it with the one measured in Et, and with the 
results of the deconfining model ||. In Section 4, we compute the J/ip suppression in pA 
and SU in the comovers approach and compare the normalization in these systems with 
the one in PbPb. The scenario for J/ip suppression resulting from the comovers model is 
described in the Conclusions and confronted with the deconfining one. 
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2 Ex - Ezdc correlation. 



In the models ||, ||, one determines the J/ip suppression as a function of the impact 
parameter, b. However, b is not measurable and the NA50 collaboration uses, as a measure 
of centrality, either Et or Ezdc- It is, therefore, necessary to determine both quantities 
as a function of b. Since Et is proportional to the multiplicity, the relation between Et 
and b results from the determination of the multiplicity at each b in the rapidity region of 
the Et calorimeter. In the comovers model || this is done || in the Dual Parton Model 
(DPM). More precisely, we put : 

Mb) = \ q N c y ° c M • (2) 

Here Ny° a is the charged multiplicity in the rapidity region of the Et calorimeter. The 
factor 1/2 is introduced because only the energy of neutrals is measured by the calorimeter. 
Thus the coefficient q is close to the average energy per particle. However, the difference 
between multiplicities of positive negatives and neutrals as well as the efficiency of the 
Et calorimeter do affect the value of q. This value can be determined from the position 
of the "knee" of the Et distribution of MB events measured by the NA50 collaboration. 
We obtain q = 0.62 GeV § (see below). 

The energy of the zero degree calorimeter is given by 



E ZDC {b) = [A-n A {b)]E in + an A {b) E in . (3) 

Here A — n A (b) is the number of spectator nucleons of A and E in = 158 GeV is the beam 
energy. While the first term in the r.h.s. of Eq. (|3|) gives the bulk of Ezdc-, the latter 
corresponds to the contamination by secondaries emitted very forward [ID| - assumed to 
be proportional to the number of participants, Here also the value of a can be 

precisely determined from the position of the "knee" of the Ezdc distribution of the MB 



event sample measured by NA50 ||10|| . We obtain a = 0.076. 

Eqs. (H) and (H) give the relation between b and E T and b and E Z dc-, respectively. 
These relations refer to average values and do not contain any information about the tails 
of the Et or E Z dc distributions. Eqs. (|]) and fl3|) also lead to a correlation between 
(average values of) E T and E Z dc- This correlation is shown in Fig. 2, and gives a good 
description of the experimental one ||. We see from Fig. 2 that the Et — Ezdc correlation 
is close to a straight lineQ and therefore can be accurately extrapolated beyond the knee 
tThis is due to the fact that Ny° a (b) in Eq. (||) is practically proportional to nA{b) (see Fig. 1 of [|J). 
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of the Et and Ezdc distributions. It turns out that this extrapolation describes the data 
quite well|. Let us discuss this point in detail. It is well known that the correlation Et — b 
can be described by a Gaussian : 



P(E T ,b) = 1 exp 

'2TiqaE T (b) 



[Et - E T {b)f 



(4) 



2qaE T (b) 

with q = 0.62 GeV and a = 0.60 ||. The resulting MB distribution is indistinguishable 



is 



from the solid line of (1998 data). The MB distribution of the 1996 data |0 
reproduced with q = 0.62 GeV and a = 0.852 §0. 

Applying the Et — Ezdc correlation resulting from Eqs. (0) and we obtain the 
Ezdc distribution of MB events shown in Figs. 3. We see that the NA50 data are well 
described, not only up to the knee, but also in the tail of the distribution. The turn- 
over in the data at large Ezdc is due to the fact that they have not been corrected for 
efficiency The comparison of Figs. 3a and 3b is quite instructive. In Fig. 3a it seems 
that the calculated curve has too much curvature and overshoots the data at small Ezdc- 
However, this is no longer the case in Fig. 3b. Here, on the other hand, there is some 
discrepancy in the region 6 < Ezdc < 12 TeV - which is not present in Fig. 3a. Finally, 
the calculated curve is slightly too broad at the tail in Fig. 3a and too narrow in Fig. 3b. 
From this comparison, we conclude that Eq. (4), supplemented with the Et — Ezdc 
correlation of Fig. 2, gives a good description of the Ezdc distribution when both the 
1996 and the 1998 data are considered. 

The consequences of this result are quite interesting. Indeed, the E T — E Z dc is essen- 
tially a correlation between multiplicity (which in the rapidity region of the calorimeter is 
practically proportional to the number of participants) and number of spectators. Hence 
it cannot be affected by the dimuon trigger[]. Therefore, not only the MB, but also the 
Ezdc distributions of J/ip and DY event samples can be obtained from the corresponding 
ones versus Et applying the Et — Ezdc correlation. This is an important result since it 
allows to relate to each other the J/ip suppression versus E T and versus E Z dc (obtained 

■''One can understand the physical origin of this extrapolation if one assumes that a fluctuation in Et 
is essentially due to a fluctuation in ua - which, in turn, produces a corresponding fluctuation in Ezdc, 
viaEq. (|). 

§ At first sight these sets of values look very different from the ones used by the NA50 Collaboration. 
Nevertheless, they reproduce the same Et distribution. This is due to the fact that the product qa, 
which, according to eq. (|]), determines the width of the distribution, is very similar in the two cases. As 
for the difference in the values of q it is just due to its definition, which is different in the two approaches 
(eq. (g), in our case). 

^Except in the tails, where the Et loss affects the Et distributions of J/tp and DY without changing 
the Ezdc ones. 
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with a different calorimeter) and check their consistency. 



3 J/ip suppression versus Ezdc- 

We show in Fig. 4 the results for the ratio of J/ip to DY cross-sections versus Ezdc 
obtained in the deconfining model and in the comovers model ||. The curves, are 
obtained from the corresponding ones versus Et (upper solid and dashed cuerves in Fig. 
1), applying the E T — E Z dc correlation given by Eqs. (f2|) and (^) and shown in Fig. 2. 
We keep the absolute normalizations unchanged. No Et loss is incorporated in these 
calculations (see below). 

Comparing the data with the model predictions, we see that a better description of 
the small Ezdc region is obtained when the Et fluctuations are taken into account. This 
was to be expected since the fluctuations in Et and Ezdc are related to each other via 
the Et — Ezdc correlation - as shown in Fig. 2. 

As discussed above, the effect of the E T loss is not included here since E Z dc is n °t af- 
fected by the dimuon triggerfj]. Nevertheless, the curve obtained in the comovers approach 
does describe the most central data while it failed to describe the data versus Et at large 
Et (see Fig. 1). This result lends support to the interpretation in ref. - namely that 
this discrepancy is due to the -E^-loss induced by the J/ip trigger. 

We also see in Fig. 4 that the fall-off for central events is stronger in the deconfining 
than in the comovers model. The present data do not allow to discriminate between them. 
However, this should be possible when the shape and absolute normalization of the data 
are better known - possibly with the 2000 NA50 data. 

Turning to the large Ezdc region (peripheral events), we see that the data favor the 
comovers model. Indeed, it is clear that, in the present data, the ratio Rmb(Ezdc) 
at large Ezdc falls more steeply than the NA50 nuclear absorption curve - fitting pA 
and SU data. If confirmed, this would clearly disfavor the deconfining model - where 
nuclear absorption is assumed to be the only source of suppression below the deconfining 
threshold. On the contrary, such a feature is expected in a comovers approach, and is 
clearly seen in Fig. 4. 

An important observation from Fig. 4 is the fact that the pattern of J/ip suppression 

"Actually, some effect of the Et loss could be present in the J/ij) suppression versus Ezdc via the 
second term of Eq. (3). However, this effect should be rather small due to the large rapidity separation 
of the dimuon trigger and the zero degree calorimeter. 
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in the NA50 data is different in Et and Ezdc ~ at least in the region 9 <Ezdc ^ 14 
TeV. This is at variance with the results in Section 2, according to which the shapes of 
the two curves should be the same. More precisely, if the shape of the J/ip suppression 
in Ezdc is different from the one in Et for whatever reason (for instance, due to the 
different calorimeters), the shape of the MB distribution in the two variables should also 
be different - and it is not. This point needs to be clarified. It should be noted here 
that the ratio of J/ip to DY differential cross-sections obtained in the standard analysis 
[ PHI , although consistent with the one obtained in the MB analysis, do not show any 
sign of a different pattern with respect to the one versus Et- This is illustrated in Fig. 5 
where the standard analysis data are compared with the same theoretical curves of Fig. 4. 
Note also that in the MB analysis the absolute normalization is not measured. It is 
determined from the one obtained with the standard analysis. This determination can 
not be precise due to the uncertainty in the shape of the J j ip suppression discussed above. 



Actually, the adjustment of the absolute normalization has been done |10| in the region 
8 < E Z cd < 17 TeV - where the difference in the shape of the E T and E Z dc distributions 
is the largest. 

Further insight in the comparison of comovers and deconfining models can be gained 
by studying the J/ip suppression in lighter systems. 



4 J /if; suppression in pA and SU. 

Let us compute the ratio J/ip over DY in SU at 200 GeV/c per nucleon in the model 
|J. We use, of course, the same values of the parameters as in PbPb : a abs = 4.5 mb 
and a co = 1 mb. To get this ratio versus b, the only new ingredient is the multiplicity of 
comovers - which is again computed in DPM, in the way described in ||. In this case, 
only R(E T ) is available (R(E Z dc) has not been measured in SU). In SU, data do not 
extend beyond the knee of the ^-distribution. Therefore, effects such as E T fluctuations 
or Et loss are not relevant here. To compute R(Et), we also need the Et — b correlation 
which is parametrized as in Eq. (|J). The parameters q and a have been obtained from a 



fit of the ^-distribution of DY given in [pT |. We obtain q = 0.69 GeV and a = 1.6. 



Our results are shown in Fig. 6a. We see that the Et dependence of the suppression 
is reproduced - in spite of the fact that our a a b s is smaller than in the NA50 nuclear 
absorption model. This is due to the suppression by comovers. Actually, the effect of the 
comovers in SU is comparatively small for all values of Et- However, its contribution 
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increases with Et and compensates for the difference in the values of a a b s in the two 
approaches. 

A recent reanalysis of the pA data at 450 GeV |L2| leads to a value a a bs = 4.7 ±0.8 mb 
in good agreement with the one used in the comovers approach. This value is significantly 
lower that the one, a a bs = 6.4 mb, used in the NA50 nuclear absorption curve. A value 
(tabs = 4 mb is used in the deconfining model of ref. ||13|| . However, within the deconfining 
approach such a value of a a b s has the drawback of producing a J/ip suppression pattern in 
SU which is flatter than the data - unless, of course, one reduces substantially the value 
of the threshold, thereby introducing some anomalous suppression in SU. 

Let us now discuss the absolute normalization of the curve in Fig. 6a which is 46.8. 
This number has to be compared with the ratio of J / if) to DY pp cross-sections, rescaled 
at 200 GeV, which is 46.6 ± 5 ||14|| . This nice agreement between pp and SU tends to 
indicate that the pA data will also be reproduced (see also [[HI). This is indeed the case, 
as shown in Fig. 6b. 

Our result gives an A-dependence between pp and pU A a with a = 0.943. This is 
to be compared with the NA38 value a = 0.919 ± 0.015 |14[ and with the E866 [16] one 
a = 0.955 ± 0.02 ± 1 % systematics. 

We can now compare the absolute normalization of our curve for pp, pA and SU (46.8) 
with the corresponding one for Pb Pb (59.4). As we see, there is a 27 % difference. Accord- 



ing to the NA50 estimates |T7| , 9 % of this difference is due to the rescaling from SU at 200 
GeV to Pb Pb at 158 GeV. This factor 1.09 contains both energy and isospin corrections]^ 
Therefore, it remains a real discrepancy of about 17 % in the relative normalization of pp, 
pA and SU, on the one hand, and Pb Pb on the other hand. This point needs clarification. 



5 Conclusions 

The deconfining - NA50 scenario, has been described in the first lines of this work. The 

alternative scenario obtained in the comovers model || is summarized in Fig. 7. The 

curves in Fig. 7 are obtained from the ones in Figs. 1 (PbPb), 6a (SU) and 6b (pp and 

**N. Armesto (private communication) has recalculated these corrections. He has confirmed the NA50 
results concerning isospin. However, he finds an energy dependence of the DY practically identical to 
the one of J/ip estimated by NA50 - leading to a ratio of J/ip to DY cross-sections practically energy 
independent. This would solve the problem with the rescaling in energy used by NA50 in which the 
central value for the ratio J/ip over DY in pp decreases between 158 GeV (48.9) and 200 GeV (46.6) and 
increases from 200 to 450 GeV (54.7). 
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pA) , using the relation between A and L (in pA) and and L (in SU and PbPb) given 
by NA38-NA50. In this Figure we see that there are two different curves for pA (dashed 
line) and SU (full line). This is due to the effect of the comovers. Contrary to nuclear 
absorption, which is a universal function of L, the comovers contribution at a given value 
of L is different in different systems. In particular, its effect turns out to be negligible in 
pA. As already seen in Fig. 6a and 6b the agreement with the data is good, both in pA and 
SU. The PbPb data are not shown in this Figure. They follow quite well the shape of the 
theoretical curve (dashed line), as demonstrated in Fig. 1 (lower dashed line). However, 
their absolute normalization is 17 % higher than the theoretical one, as discussed above. 
In the comovers model, the J/ip suppression in PbPb is larger than the one in SU (at 
the same L) and significantly steeper - indicating an important anomalous suppression, 
present even in the most peripheral events. The data do support the latter result - which 
favors the comovers scenario. For central events, it lends support to the interpretation of 
ref. ||, according to which the discrepancy between the comovers model and Rmb, eq. 
(HD, is due to the E T loss induced by the J ftp trigger. Indeed, the data versus E Z dc (which 
are not affected by this Et loss) are well reproduced. Hopefully, a clear-cut discrimination 
between the deconfining and comovers scenarios will be possible with the 2000 NA50 data. 
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Figure Captions 



Fig. 1 : Ratio of J/ip to DY cross-sections versus E? in Pb Pb collisions at 158 GeV 
per nucleon in the deconfining model [|J (solid) and in the comovers model || (dashed). 
In both cases the upper (lower) curves are obtained without (with) Et loss. 

Fig. 2 : The correlation Et — Ezdc obtained from Eqs. (|2|) and (|3|) and its extrapolation 
to the tails of these distributions (see main text). 

Fig. 3 : E Z dc distribution of MB events obtained from Eq. (^) applying the E t — E Z dc 
correlation of Fig. 2. The data in Fig. 3a were taken in 1998 [JTo|] and those in Fig. 3b 
were taken in 1996 [HJ- 



Fig. 4 : Ratio of J/ip to DY cross-sections versus Ezdc m Pb Pb collisions at 158 GeV 
per nucleon in the deconfining model [3J (solid) and in the comovers model || (dashed). 
These curves are obtained from the corresponding ones versus Ex (upper solid and dashed 
lines in Fig. 1) applying the E T — E Z dc correlation of Fig. 2. The dotted line is obtained 
without Et fluctuations. The data, obtained with the minimum bias analysis, are from 



ref. JII|. The NA50 nuclear absorption curve [ID| is also shown. 



Fig. 5 : Same theoretical curves as in Fig. 4, compared with the data obtained in the 
standard analysis |T(J . 



Fig. 6a : The ratio of J/if> to DY cross-sections as a function of Et in SU collisions at 
200 GeV per nucleon obtained in the comovers model ||. The data are from |TT |. 



6b : The ratio J / if) over A in pp and pA collisions at 450 GeV as a function of A obtained 
in the comovers model ||. The data are from [14 . 



Fig. 7 : The ratio of J/ip to DY cross-section versus L || for pp, pA (dotted line) 
and SU at 200 GeV (solid line) and PbPb (dashed line) at 158 GeV, obtained in the 
comovers model [El with a common normalization 46.8. Note that this normalization is 



not calculable in the model. The data are from (TT| [|1J . 
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Fig. 1 
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Fig. 2 
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Fig. 3a 
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Fig. 3b 
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Fig. 4 
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Fig. 5 
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Fig. 6 
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Fig. 7 
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